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The existence of supermassive black holes is supported by 
a growing body of observations. Supermassive black holes 
and their formation events are likely candidates for detection 
by proposed long- wavelength, space-based gravitational wave 
interferometers like LISA. However, the nature of the progen- 
itors of supermassive black holes is rather uncertain. Super- 
massive black hole formation scenarios that involve either the 
stellar dynamical evolution of dense clusters or the hydrody- 
namical evolution of supermassive stars have been proposed. 
Each of these formation scenarios is reviewed and the evolu- 
tion of supermassive stars is then examined in some detail. 
Supermassive stars that rotate uniformly during their secu- 
lar cooling phase will spin up to the mass-shedding limit and 
eventually contract to the point of relativistic collapse. Super- 
massive stars that rotate differentially as they cool will likely 
encounter the dynamical bar mode instability prior to the 
onset of relativistic collapse. A supermassive star that under- 
goes this bar distortion, prior to or during collapse, may be a 
strong source of quasiperiodic, long-wavelength gravitational 
radiation. 



I. SUPERMASSIVE BLACK HOLES 

There is a large body of observational evidence that 
supermassive black holes (SMBHs) exist in the centers 
of many, if not most galaxies (see e.g., the reviews of 
Rees fij and Macchetto 0). The masses of SMBHs in 
the centers of more than 45 galaxies have been estimated 
from observations Q and there are more than 30 galaxies 
in which the presence of a SMBH has been confirmed M . 
The properties of a few confirmed SMBHs are given in 
Table 1. 

Several properties have been deduced from these ob- 
servations |Q,||J|-|l|]. A correlation between the SMBH 
mass and both the mass and the velocity dispersion of the 
bulge of the host galaxy is observed. These results sug- 
gest that the formation and evolution of SMBHs and the 
bulge component of their host galaxies may be closely re- 
lated. The largest SMBHs arc found in elliptical galaxies 
and thus may result from galaxy mergers, as ellipticals 
themselves are thought to form via merger. 

Because of their size and mass, SMBHs are expected 
to be sources of long- wavelength, low-frequency gravi- 
tational radiation. Thus, phenomena involving SMBHs 
might be detectable with proposed space-based gravita- 
tional wave detectors, like the Laser Interferometer Space 



Antenna (LISA) [nsl. For example, LISA might be able 
to detect the collapse of a supermassive star to a SMBH, 
the coalescence of two SMBHs, or the coalescence of a 
compact star and a SMBH [n#-n8|. The event rates of 
these phenomena are uncertain due to the uncertainty in 
the nature of the progenitors of SMBHs. However, the 
rates could be appreciable because SMBHs are present in 
many galaxies. 



A. Formation Mechanisms 

The possible formation mechanisms for SMBHs involve 
either the stellar dynamics of dense star clusters or the 
hydrodynamics of supermassive stars (SMSs). 



1. Stellar Dynamical Routes 

In general terms, dense star clusters are formed via a 
conglomeration of stars, produced by fragmentation of 
the primordial gas. In one cluster scenario for SMBH 
formation, massive stars form via stellar collisions and 
mergers in the cluster and then evolve into stellar-mass 
black holes. The merger of these holes, as they grow and 
settle to the cluster center, leads to the build up of one 
or more SMBHs |l|,|2^,|l|l . 

An alternative cluster scenario centers on the catas- 
trophic collapse of the cluster core. Shapiro and Teukol- 
sky have numerically demonstrated that cluster cores can 
become dynamically unstable to relativistic collapse [£1| . 
They also showed that if the core is made up of com- 
pact stars, it will evolve to the onset of collapse on a 
relaxation timescale [p2|. As the compact core evolves 
to this point of relativistic instability, the mass of its 
black holes will increase through collisions and gravita- 
tional radiation driven mergers. This will accelerate the 
approach to relativistic instability, leading to SMBH for- 
mation [PII2I. 



2. Hydrodynamical Routes 

SMSs may contract directly out of the primordial gas, 
if radiation and /or magnetic field pressure prevent frag- 
mentation p5|-|27|] (see also p^-|30(|). Alternatively, they 



may build up from fragments of stellar collisions in clus- 
ters gH. 

The evolution of SMSs will ultimately lead to the onset 
of relativistic instability pJ^-jSq] . If the mass of the star 
exceeds lO^M©, the star will then collapse and possibly 
form a SMBH. If the star is less massive, nuclear reactions 
may lead to explosion instead of collapse. 



II. SUPERMASSIVE STARS 



The major result of Baumgarte and Shapiro's work 
is that the following universal ratios exist for the criti- 
cal configuration at the onset of relativistic instability: 
T/\W\, R/M, and J/W^. Here J is the total angular 
momentum of the star. These ratios are completely in- 
dependent of the mass of the star or its prior evolution. 
Because uniformly rotating SMSs will begin to collapse 
from a universal configuration, the subsequent collapse 
and the resulting gravitational waveform will be unique. 



Identifying the scenario by which SMBHs form is of 
fundamental importance to a number of areas of astro- 
physics. The remainder of this paper examines one of the 
possible progenitors of SMBHs discussed above, SMSs. 

Supermassive stars are radiation dominated, isen- 
tropic, and convective |^,^,g8[ and are thus well repre- 
sented by a n = 3 polytrope. As mentioned previously, 
if the star's mass exceeds lO^M©, nuclear burning and 
electron/positron annihilation are not important. 

After formation, a SMS will evolve through a phase of 
quasistationary cooling and contraction. If the SMS is ro- 
tating when it forms, conservation of angular momentum 
requires that it spin up as it contracts. The evolutionary 
path taken by the SMS during this cooling phase depends 
on the strength of its viscosity and magnetic fields and 
on the nature of its angular momentum distribution. 



A. Cooling Evolution of Uniformly Rotating 
Supermassive Stars 

There are two possible evolutionary regimes for a cool- 
ing SMS. In the first regime, viscosity or magnetic fields 
are strong enough to enforce uniform rotation through- 
out the star as it contracts. 

Baumgarte and Shapiro Q have recently studied the 
evolution of a uniformly rotating SMS up to the onset 
of relativistic instability. They demonstrated that a uni- 
formly rotating, cooling SMS will eventually spin up to 
its mass shedding limit. The mass shedding limit is en- 
countered when matter at the star's equator rotates with 
the Keplerian velocity. The limit can be represented as a 
ratio Pshed = [T /\W\)shed of the star's rotational kinetic 
energy T to gravitational potential energy W . In this 
case, j3shed = 9 x 10~^. The star will then evolve along 
a mass shedding sequence, losing both mass and angular 
momentum. It will eventually contract to the point of 
relativistic instability. 

Baumgarte and Shapiro used both a 2nd order, post- 
Newtonian approximation and a fully general relativistic 
numerical code to determine that the onset of the insta- 
bility occurs at a ratio of R/M ^ 450, where R is the 
star's radius. Here and throughout this paper G = c — 1. 
Note that a 2nd order, post-Newtonian approximation 
was needed because rotation stabilizes the destabilizing 
role of nonlinear gravity at the 1st post-Newtonian level. 



1. Collapse Outcome and Gravitational Radiation Emission 

The outcome of SMS collapse can only be determined 
with numerical, fully relativistic three-dimensional hy- 
drodynamics simulations. To date, such simulations have 
only been published for nearly spherical collapse. The nu- 
merical simulations of Shapiro and Teukolsky p9| demon- 
strate that this type of collapse is nearly homologous. 

In this case the collapse time TcoU is roughly the free- 
fall time at the horizon (where R = 2M) 
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in the middle of LISA's frequency band of 10 — 1 Hz 

The amplitude h of this burst signal can be estimated 
roughly in terms of the star's quadrupole moment 



h<e- 



2M2 



< e5 X 10"^^ [M/lO^Af0][d/lGpc]- 
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Here d is the distance to the star and e ~ r/|iy| is a 
measure of the star's deviation from spherical symmetry. 
In this case, e will be much less than one even near the 
horizon, since the collapse is nearly spherical. 

Three-dimensional, general relativistic simulations of 
nonspherical, rotating collapse are underway |^^. The 
results are not available as yet. However, there are two 
possible outcomes of this type of collapse that can be 
discussed. 

The first outcome is direct collapse to a SMBH, from 
the onset of instability. In this case e will be on the 
order of one near the horizon. Thus the peak ampli- 
tude (see equation 2) of the burst signal will he h ^ 
5 X 10-1^ [M/lO6M0][d/lGpc]"i. 

Alternatively, the star may encounter the dynamical 
bar mode instability. The bar mode is the strongest of 
a set of global nonaxisymmetric instabilities that may 
be encountered by a rapidly rotating object. This in- 
stability will deform the star into a bar-shaped config- 
uration, making it a strong source of long- wavelength. 



quasiperiodic gravitational radiation. During the devel- 
opment of the instability, angular momentum and mass 
will be transported outwards ||4l|-^. This could hasten 
the star's eventual collapse. Previous linear and nonlin- 
ear analyses indicate that f3bar ^ 0.27 is likely to be an 
upper limit for the onset of this instability, for a wide 
range of polytropic equations of state and rotation laws 
(see, e.g., [|4|||). 

Baumgarte and Shapiro Bq] have estimated that 
a uniformly rotating SMS will reach /3~0.27 when 
R/M=15. This estimate awaits confirmation from three- 
dimensional, general relativistic hydrodynamical simula- 
tions. 

The frequency of the quasiperiodic gravitational radi- 
ation emitted by the bar can be estimated in terms of its 
rotation frequency to be 



few = 2fi 



fGM\ 
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2 X IQ-^RzIM/IO^Mq]-^, 
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when R/M=\b. The corresponding amplitude of the ra- 
diation, again estimated in terms of the star's quadrupole 
moment, is 



2M2 

6 X 10-1* [M/10SA'fQ][d/lGpc]-^ 



(4) 



The emission of gravitational radiation will continue as 
long as the collapsing SMS's bar shape persists. 



B. Cooling Evolution of Differentially Rotating 

Supermassive Stars 

In the opposite evolutionary regime, neither viscosity 
nor magnetic fields are strong enough to enforce uniform 
rotation throughout the cooling SMS as it contracts. In 
this case it has been shown that the angular momen- 
tum distribution is conserved on cylinders during con- 
traction iQ. Because viscosity and magnetic fields are 
weak, there is no means of redistributing angular momen- 
tum in the star. So even if the star starts out rotating 
uniformly, it cannot remain so. 

The star will then rotate differentially as it cools and 
contracts. In this case, the subsequent evolution de- 
pends on the star's initial angular momentum distribu- 
tion, which is largely unknown. 

One possible outcome is that the star will spin up to 
mass shedding (at a different value of (3shed than a 
uniformly rotating star) and then follow an evolutionary 
path that may be similar to that described by Baumgarte 
and Shapho ^. 

The alternative outcome is that the star will encounter 
the dynamical bar instability prior to reaching the mass 
shedding limit. As discussed above, the bar shape in- 
duced by this instability will make the star a source 



of quasiperiodic, long-wavelength gravitational radiation. 
The outward transport of angular momentum that occurs 
during the instability could hasten the eventual collapse 
of the star. 

New and Shapiro B^ have investigated the evolution 
of differentially rotating SMSs. Because the angular mo- 
mentum distribution of these stars is unknown, they ex- 
amined SMS models with several different initial angular 
momentum distributions. The goal of their work was to 
determine if the bar instability or mass shedding limits 
are reached by these SMS models. 

New and Shapiro's strategy was to examine equilib- 
rium sequences of SMS models, each of which was con- 
structed with a different rotation law. The individual 
models on each sequence were constrained to have the 
same M and J, since these quantities are conserved dur- 
ing the cooling evolution of a SMS. However, the models 
along a sequence have decreasing entropy and thus de- 
creasing axis ratio Rp/ R^, where Rp is the polar radius 
and Re is the equatorial radius. A sequence is thus rep- 
resentative of the quasistatic, cooling/contracting evolu- 
tion of a single SMS. This quasistatic approximation is 
appropriate because the cooling timescale is much longer 
than the hydrodynamic timescale for M <^ IQ^^Mq [Esj. 
New and Shapiro examined each of the sequences to de- 
termine whether the limits P^ar or Psted are reached for 
a SMS model with the given rotation law. 

All of the sequences New and Shapiro examined were 
constructed with Hachisu's self-consistent field (HSCF) 
technique Q. The HSCF method builds individual 
models in hydrostatic equilibrium, such that their pres- 
sure, gravitational, and centrifugal forces are in balance. 

The HSCF method requires the choice of a barotropic 
equation of state P — P{p). All of the SMS models 
examined by New and Shapiro were constructed with a 
polytropic equation of state for which 



P^ Kp 



l+l/n 



(5) 



As mentioned above, the structure of a SMS is well repre- 
sented by an 71 = 3 polytrope. In this equation of state, 
the polytropic constant K is a measure of the specific 
entropy of the model. 

The selection of appropriate rotation laws was con- 
strained by the fact that each rotation law chosen had 
to enforce conservation of the contracting star's specific 
angular momentum profile |I6[ . 

One rotation law that satisfies this constraint is the 
so-called n'=3 law Efl. Its specific angular momentum 
profile j{m) is 



j{m) =01+02 1 - 



M 



03 1- 



m{w) 
M 



(6) 



where M is the total mass of the system, m is the mass 
interior to cylindrical radius w, and the numerically de- 
termined constants are oi=13.27, 02=163.3, 03=-176.5, 



j{m) = (1 + q){J/M) [1 - (1 - mivj)/My^^] 



a2=0.2353, and a3=0.2222. This angular momentum 
profile is identical to that of a uniformly rotating, spher- 
ical n—3 polytrope. Thus an equilibrium sequence con- 
structed with this rotation law is representative of the 
evolution of a SMS that rotates uniformly prior to its 
cooling phase. New and Shapiro [^ have constructed 
an equilibrium sequence with this rotation law. 

Because the initial rotation profiles of SMSs are un- 
known, New and Shapiro also examined sequences con- 
structed by Hachisu, TohHne, and Eriguchi ( |^; here- 
after, HTE). These sequences were built with the follow- 
ing parameterized angular momentum profile: 

(7) 

Here, the index q specifies the rotation law. Note that 
the limiting case of q=0 corresponds to the j-constant 
rotation law. 

Nearly spherical models built with these "q-laws" are 
differentially rotating. Thus, HTE's n=3, q-law se- 
quences are representative of the evolution of SMSs, with 
a wide range of initial differential rotation profiles. 



1. Evolutionary Scenarios 

A detailed discussion of the properties of the n'=3 se- 
quence, constructed by New and Shapiro, and the q-law 
sequences of HTE can be found in |^^ . In what follows 
we summarize these properties as they relate to the evo- 
lutionary scenarios of differentially rotating SMSs. 

Density contour plots of selected models from the n'—i 
sequence are shown in Figure 1. The n'=3 sequence ter- 
minates due to mass shedding at Pshed ^ 0.30. Pshed 
exceeds the likely upper limit for the bar instability 
fibar <i 0.27. Thus we expect that a SMS with this rota- 
tion law should never reach the mass shedding limit, but 
will instead encounter the dynamical bar instability near 
Rp/Re ~ 0.004. This is the axis ratio of the model with 
/3 = 0.27 (see Figure Id). 

No mass shedding limits exist on HTE's q-law se- 
quences. Each of these sequences makes a continuous 
transition from spheroidal to toroidal configurations at 
values of Ptrans > 0.33 > 0.27 ^ Pbar- Thus, n=3 mod- 
els with these q-indexed laws would likely encounter the 
bar mode as spheroids. 

We note that the hydrodynamical study of B3] indi- 
cates that Pbar may be less than 0.27 for rotation laws 
that place a significant amount of angular momentum in 
equatorial mass elements. Their results predict that the 
m=2 stability limit may be less than ~ 0.20 for models 
with the n'—3 rotation law. Their simulations, of 7i=1.5 
polytropes, also suggest that a one-armed spiral, m—1 
mode may become increasingly dominant over the m—2 
mode as the equatorial concentration of angular momen- 
tum increases. Note that the grid resolutions used in H] 
were likely not sufficient to accurately model the develop- 
ment of instabilities in models with these extreme differ- 
ential rotation laws Q|. However, the results of Q] and 



the linear and nonlinear stability analyses of [Q confirm 
that 0.27 is likely to be an approximate upper limit to 
Pban for a variety of polytropic indices and rotation laws. 
The analysis presented in |^ assumes that the m=2 bar 
mode is the dominant mode and that (3bar <, 0.27. 

Even if the actual value of /3bar is less than 0.27, the 
qualitative nature of the results of New and Shapiro 
would not change. That is, the sequences they exam- 
ined would still have models that are unstable to the bar 
mode. In addition, their quantitative estimates of the 
characteristics of the gravitational radiation emission 
presented below would only change by a numerical factor 
of order 1-10, even if (3bar were as low as 0.22. 



2. Gravitational Radiation Emission 

The results of New and Shapiro |47) indicate that a 
bar mode phase is likely to be encountered by differen- 
tially rotating SMSs with a wide range of initial angular 
momentum distributions. Hydrodynamical simulations 
are needed to follow the evolution of the star as the in- 
stability develops, to compute the gravitational radiation 
waveforms emitted, and to determine the fate of a SMS 
that undergoes the bar instability. 

Previous hydrodynamical simulations of the bar insta- 
bility in models with other rotation laws and a different 
polytropic equation of state indicate that the outcome 
of this instability is a persistent bar-like structure that 
emits quasiperiodic gravitational radiation over many 
cycles [ 40| , p^ 43,5^]. If a similar outcome results from 
the bar instability in a SMS, the quasiperiodic, long- 
wavelength gravitational radiation emitted could be de- 
tected by LISA. 

The frequency of these quasiperiodic gravitational 
waves can be estimated from the expected bar rota- 
tion rate ^bar- The model on New and Shapiro's 
n'=3 sequence with (3 '^(3bar=0 -27 has a central rota- 
tion rate 17^=1.02 x lO^^Hz [MePztRijV^^- Here Mq = 



Af/lO^Mo, i?i7 = (i?e)o/10^^cm, and /5_5 = Po/W- 
The subscript denotes the value is for the (nearly) 
spherical progenitor star, prior to the start of its cool- 
ing phase. In previous hydrodynamics simulations of the 
bar mode instability |^^, Qbar was '^ OAQc- With this 
relation between flbar and Oc, the gravitational wave fre- 
quency few can be estimated to be 



few = 2fbar 



n 



bar 



27r 



IT 



1 X lO^^Hz [MePZiRu] 



(8) 



For a SMS of IO^Mq, which begins as a slowly rotating 
star of radius lO^^cm with /3o=10"^, this yields a fre- 
quency of 1 X 10~^Hz. This frequency is in the range in 
which LISA is expected to be most sensitive, 10~^-lHz. 
The choice /3q=2 x 10~''/i?i7 is the maximum value of 
Po for which few (=10~^Hz) is still in EISA's range of 
sensitivity. 



The strength of the gravitational wave signal can be 
estimated roughly to be 



GQ GMRlJ, 

C4 d c4 



barJ bar 



4 X 10" 



IGpc 



2o-l D-1 



^hP-s^n 



(9) 



where Q is the second time derivative of the star's 
quadrupole moment and d is the distance, which we 
scale to 1 Gpc (the Hubble distance is ~ 3 Gpc). Here 
we have used Rbar ^ -Re=4.14 x 10^^/3_5i?i7 cm for the 
n'=3 model of [^ that reaches the point /3 '^ (3bar- 

The bar will decay on a secular timescale due to dis- 
sipative effects. For differentially rotating SMSs, the 
largest source of dissipation will be gravitational radia- 
tion. The gravitational radiation damping timescale TQ\y 
is approximately 



tgw 



V dt Jgw 



(10) 



where T is the rotational kinetic energy and {dE / dt)Q\Y 
is the rate at which gravitational radiation carries energy 
away from the system. Recall that T — (3\W\. The 
gravitational potential energy \W\ ~ GM^/Ri,ar- Thus, 



T = /3 



GM'^ 
R 



The radiation rate can be estimated as 



,dE. 
y dt'GW 



G M 
l^R^' 



(11) 



(12) 



In this case the characteristic velocity of the system is 
V = (/3GM/i?)^/^. Substitution of equations (11) and 
(12) into equation (10) yields 



Tgw 



Rt 



G3 ^2^3 
lxl04yrsAfg-3[/3_5i?l7]''. 



(13) 



The number of cycles M for which the signal will per- 
sist is 



M ^ Tgw Jgw 

^4xl09[A/6/3_5i?i7] 



-5/2 



(14) 



The quasiperiodicity of such a signal will assist in its 
detection [1^. 

The fraction of the mass {dM /M)gw radiated via 
gravitational radiation over the interval tgw can be es- 
timated as 



fdM\ 
VW) 



,dE. 
^Gw[-^) 



1 



GW 



dt'GW Mc^ 



G fiM 

lP~R 

1 X 10-3Af6[/3_5i?i7]-i. 



(15) 



III. CONCLUSIONS AND FUTURE WORK 

There is strong evidence that SMBHs exist. However, 
the nature of their progenitors is uncertain. Proposed for- 
mation mechanisms involve the evolution of stellar clus- 
ters or SMSs. 

Recent studies indicate that differentially rotating 
SMSs, and possibly uniformly rotating SMSs as well, are 
likely to encounter the dynamical bar instability and thus 
emit quasiperiodic, long-wavelength gravitational radia- 
tion. SMSs will also emit burst gravitational wave signals 
as they undergo relativistic collapse. 

Linear and nonlinear stability analyses are needed to 
determine precisely the onset of the m=2 dynamical bar 
instability for the SMS models considered here. Such 
analyses are also needed to determine the relative im- 
portance of various unstable nonaxisymmetric modes in 
SMS models (as a dominant m—1 mode would change the 
characteristics of the gravitational radiation emission). 

In addition, three-dimensional hydrodynamical simu- 
lations are necessary to study the evolution of SMSs as 
they undergo the bar instability and/or relativistic col- 
lapse, to compute the gravitational waveforms emitted, 
and to determine the final fate of the star. 

Future investigations involving hydrodynamical sim- 
ulations of SMS models would benefit from improved 
knowledge of initial conditions, such as an appropriate 
value for Pq. These appropriate initial conditions could 
be determined from studies of large-scale structure and 
cosmology. 

Most importantly, studies are required to assess the 
roles of viscosity and magnetic fields in rotating SMSs to 
judge whether they are sufRcient to drive these configu- 
rations to uniform rotation prior to bar instability. 
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FIG. 1. Snapshots of a cooling SMS. Density contours of 
selected models on the n'=3 equilibrium sequence are shown 
in the (2; > 0, z > 0) plane. The maximum density is nor- 
malized to unity. The highest density contour level is 0.9; 
subsequent contour levels range from 10~^ to 10~^° and are 
separated by a decade. The axis ratios Rp/Rs of the mod- 
els displayed are (a) 1.00; (b) 0.700; (c) 0.300; (d) 0.004; 
(e) 0.002. The model with Rp/Re=0.QQ4, shown in (d) has 

f3=Pbar ~ 0.27. 



Galaxy 


Mass (Mq) 


Radius (pc) Reference 


NGC 4258 


~ 3.6 X 10' 


< 13 
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M87 


- 3.2 X 10^ 


< 3.5 [ 


1 


NGC 4151 


^10« 


<60 [ 


1 


Milky Way 


2.6 - 3.3 X 10"^ 


< 8 


J 



TABLE I. The host galaxies, masses, and radii of selected 
confirmed SMBHs. 
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